In silico and in vitro studies confirm Ondansetron as a novel acetylcholinesterase and butyrylcholinesterase inhibitor

Alzheimer’s disease (AD) is a progressive neurodegenerative disorder that is growing rapidly among the elderly population around the world. Studies show that a lack of acetylcholine and butyrylcholine due to the overexpression of enzymes Acetylcholinesterase (AChE) and Butyrylcholinesterase (BChE) may lead to reduced communication between neuron cells. As a result, seeking novel inhibitors targeting these enzymes might be vital for the future treatment of AD. Ondansetron is used to prevent nausea and vomiting caused by chemotherapy or radiation treatments and is herein shown to be a potent inhibitor of cholinesterase. Comparison is made between Ondansetron and FDA-approved cholinesterase inhibitors Rivastigmine and Tacrine. Molecular docking demonstrates that interactions between the studied ligand and aromatic residues in the peripheral region of the active site are important in binding. Molecular dynamics simulations and binding pose metadynamics show that Ondansetron is highly potent against both enzymes and far better than Rivastigmine. Inhibitor activities evaluated by in vitro studies confirm that the drug inhibits AChE and BChE by non-competitive and mixed inhibition, respectively, with IC50 values 33 µM (AChE) and 2.5 µM (BChE). Based on the findings, we propose that Ondansetron may have therapeutic applications in inhibiting cholinesterase, especially for BChE.

Ondansetron (trade name Zofran) is a type 3 serotonin (5-hydroxytryptamine) receptor (5-HT3) antagonist ( Fig. 2A), which is prescribed against vomiting and nausea caused by chemotherapy or radiation treatment 14,15 . In comparison with other antiemetic drugs, the side effects are moderate, with high safety and efficacy. It is also used as an antiemetic drug for the first trimester of pregnancy, as well as in relief of opioid-and alcohol-induced pruritus, anxiety disorders, withdrawal syndrome, and gastrointestinal motility disorders 16,17 .
Some studies have explored the effect of Ondansetron on cholinesterase inhibition. For example, the combination of Ondansetron with Pyridostigmine as an organophosphorus pretreatment compound showed significant decrease in AChE activity in red blood cells of guinea pigs 18 .
It has also been demonstrated that the 5-HT3 receptor antagonist Ondansetron together with the FDA approved drugs Donepezil potentiates the effects of AChE inhibition on the neuronal network oscillations in the rat dorsal hippocampus 19 .
Another usage of Ondansetron as a therapeutic agent is in the treatment of psychiatric disorders which involve abnormalities of interception and associated neural circuitry centered on the insula. Different doses of Ondansetron were applied for patients whereby it was finally suggested that 24 mg of this drug modulates the hyperactivity in these regions 20 .
Regarding BChE inhibition, recent in silico and in vitro studies of a new series of Tacrine derivatives (spiro[chromeno [4,3-b]thieno [3,2-e]pyridine]-7-amines) showed that these exerted good BChE inhibitory activity, thereby proving to be promising candidates for evaluation in synthetic AD models 21 .
Based on the above findings, it is clear that Ondansetron can affect the nerve system. As a result, this study aims to explore Ondansetron as an inhibitor of cholinesterases. As a first step, different in silico approaches were applied to investigate the interaction of Ondansetron with AChE and BChE. Comparison was made between this   www.nature.com/scientificreports/ drug in both systems with the two FDA-approved cholinesterase inhibitors Rivastigmine and Tacrine (Table 1) to evaluate its efficacy. In vitro and kinetics studies of the inhibition strongly confirm the computational model. The calculation of K i and IC 50 values for these enzymes along with Arrhenius plots verify the drug binding and modes of action, and indicates that Ondansetron has stronger affinity towards BChE compared to AChE. Moreover, our study shows that Ondansetron is more potent than Rivastigmine in both enzymes but less so compared to Tacrine. For the latter, however, severe adverse effects including hepatotoxicity has led to its withdrawal in many countries 11 ; thus new and safer drugs targeting cholinesterases is of significant importance. To validate the desired properties of Ondansetron, in vivo studies must be performed to confirm its strength, reliability, longevity, efficiency, and potency.

Materials and methods
In silico studies. Protein preparation. The three-dimensional crystal structures of the target enzymes AChE (PDB ID:4M0E) 22 and BChE (PDB ID:5DYW) 23 were downloaded from the Protein Data Bank (PDB) (http:// www. rcsb. org). The proteins were prepared and refined using the Protein Preparation Wizard in Maestro (Schrödinger 2021-4, http:// www. schro dinger. com). Bond orders were assigned during the preprocessing stage of the crystal structures, and after retrieving missing loops or side chains, all water molecules beyond 3.0 Å were deleted from the system. Protein hydrogen bond assignments were optimized and protonation states at pH 7 determined using PROPKA 24 . Finally, a restrained minimization with the OPLS4 force field 25 was performed with an RMSD convergence of heavy atoms of 3.0 Å.
Ligand preparation and docking. All ligand structures were downloaded from the PubChem database (http:// www. pubch em. ncbi. nlm. nih. gov) and transferred into Maestro Schrödinger using LigPrep. The Epik module 26 was used for possible ionization states at physiological pH 7.0 ± 2.0, and the OPLS4 force field was selected for the optimization. In order to perform the molecular docking and predict the interactions of the protein-ligand complexes, the Schrödinger Induced Fit Docking (IFD) methodology was used 27,28 . All three inhibitors were docked towards the active sites of both AChE and BChE. The grid box for the docking was defined at the centroid of the binding site. During the initial docking procedure, the van der Waals scaling factor was set at 0.5 for both receptor and ligand. The Prime refinement step was set on side chains of residues within 5 Å of the ligand. No constraints were applied, and all remaining parameters were set to default.
PASS prediction. The PASS (Prediction of Activity Spectra for Substances) predictions of Ondansetron, Rivastigmine and Tacrine were carried out using the PASS-Way2Drug server (http:// www. pharm aexpe rt. ru/ passo nline/), using canonical SMILES from the PubChem server (https:// pubch em. ncbi. nlm. nih. gov/) 29 . In the PASS predictions, the probability to be active (Pa) was kept greater than 0.3, and compared with the probability to be inactive (Pi). Four possible biological activities were predicted.

Molecular dynamics (MD) simulations.
The stabilities of the ligand-protein complexes were investigated through MD simulations for 200 ns, using the Desmond engine 30 in Schrödinger (Schrödinger 2021-4, http:// www. schro dinger. com). The TIP3P force field was used to model water molecules 31 . Periodic boundary conditions were applied with a 10 Å water buffer around the protein in a cubic simulation box. Na + and Cl − ions were added to neutralize the system and to give a final NaCl concentration of 150 mM. The OPLS4 force field was used for the proteins and ligands. The isothermal-isobaric (NPT) ensemble was used, and for adjusting the temperature and pressure of the systems, the Nose Hoover thermostat 32 and the Martyna-Tobias-Klein barostat 33 were employed at 300 K and 1.01325 bar, respectively. All data analyses such as calculation of root mean square deviations (RMSD), root mean square fluctuation (RMSF), and protein-ligand contacts were obtained from the simulation interaction diagram (SID) program in Schrödinger 2021-4 (http:// www. schro dinger. com).
Clustering. To obtain representative structures for the Binding Pose Metadynamics simulations (BPMD), Desmond Trajectory clustering in Maestro (Schrödinger 2021-4) was used, based on the obtained MD trajectories. The lowest energy structures from the most populated clusters were selected for the BPMD calculations.
Binding pose metadynamics (BPMD). Binding pose metadynamics was performed using a set biasing force, to explore how stable ligands are in the binding pocket of the receptor. Weaker ligands will experience higher fluctuations with larger RMSD values in comparison with the more stably bound ones 34 . BPMD was implemented using the Desmond engine 30 in Schrödinger 2021-4 (http:// www. schro dinger. com), with the OPLS4 force field. Before the actual metadynamics run, the system was solvated in a box of TIP3P water molecules followed by several minimization and restrained MD steps to allow the system to slowly reach the desired temperature of 300 K, as well as releasing any bad contacts and/or strain in the initial starting structure 33 . 10 independent BPMD simulations of 10 ns were performed for each system, using the Root-Mean-Square Deviation (RMSD) of the ligand heavy atoms relative to their starting position as the Collective Variable (CV).
Experimental studies. Sample collection. A blood sample was obtained from a healthy volunteer. The protocol was approved by the Human Ethics Committee of Shahid Beheshti University (SBU) of Iran with the approval number of IR.SBU.REC.1401.050. All research was performed in accordance with relevant guidelines/ regulations, and informed consent was obtained from all participants. The volunteer had no significant medical disorder, was not taking any previous medication for at least  www.nature.com/scientificreports/ abuse, and no recurrent or a past history of psychiatric illness. From the volunteer, 5 ml of blood was collected in vacutainer tubes.
Acetylcholinesterase preparation. AChE of the erythrocytes is bound to the plasma membrane. The blood sample was treated with EDTA (1.5 mg/ml) as anticoagulant agent and centrifuged at 3000g for 5 min to precipitate the erythrocytes. The plasma and buffy coats were removed and discarded. Erythrocytes were washed twice with isotonic solution (NaCl 0.15 M), each time centrifuged to precipitate the intact cells. Erythrocytes were lysed by adding a hypotonic solution (NaCl 0.01 M). The solution was centrifuged at 15,000g for 20 min and the supernatant, which contained hemoglobin and other cell contents, was removed. The precipitate was washed with phosphate buffer 0.1 M, pH 7, and centrifuged at 15,000g for 20 min. The final precipitate (cell membrane) was collected for the AChE assay.
Butyrylcholinesterase preparation. BChE is present in the blood serum. The blood sample was allowed to clot without adding any anticoagulants, and then centrifuged at 3000g for 5 min. The serum as a source of BChE was separated and stored at − 20 °C for further use.
Enzyme assay. Both AChE and BChE were assayed according to the Ellman colorimetric method 35 . In brief, the assay tube contained 1500 µl 0.1 M phosphate buffer at pH 7, 50 µl 50 mM acetylthiocholine, 50 µl 10 mM DTNB (5,5′-dithiobis-2-nitrobenzoic acid), and 20 µl of either blood serum as a source of BChE, or plasma membrane solution as a source of AChE. The assay was performed either in the absence or presence of Ondansetron in the concentration range 8-84 µM. The yellow reaction product TNB (5-thio-2-nitrobenzoic acid) was monitored at 410 nm using a UV-Visible 1240 Shimadzu spectrophotometer and the activities of the enzymes calculated using the TNB extinction coefficient 13,600 M −1 cm −1 . The effect of pH and different temperatures on the enzyme activity was also measured in the absence or presence of the drug with the final concentration of 21 µM. Lineweaver-Burk plots were used to determine the kinetic parameters and inhibition types. Arrhenius plots were used to compare the activation energies of the enzymes in the presence or absence of the drug by calculating Vmax of the enzymes at different temperatures 36,37 .
The protein content of the samples was measured using the Lowry method 38 , with casein used for the standard curves.

Results and discussion
In silico studies. Molecular docking. The structure of AChE contains several important parts such as the esteratic site which harbors the catalytic triad Ser203, His447, and Glu334 (Fig. 2B), located in a 20 Å deep narrow gorge containing several conserved amino acids. A second part is the peripheral anionic site (PAS) which extends beyond Tyr337 at the catalytic/peripheral site interface to the entrance of the gorge, and contains several aromatic side chains (e.g., Tyr72, Trp86, Tyr124, Phe295, Tyr337 and Phe338). Kinetic and thermodynamic studies have shown that inhibitors can interact with either or both of the two binding regions 22 . In BChE, the catalytic triad located at the bottom of the 20 Å gorge is made up by Ser198, His438, and Glu325 (Fig. 2C) 39 . Four aromatic residues in the peripheral site, Trp82, Trp231, Tyr 332 and Phe329, plus Asp70, have been found to be important for ligand binding 40,41 .
Alignment was made between the amino acid sequences (Fig. S1), aiming to compare the active site residues in both enzymes using BChE as reference. The sequence identity and similarity between the two proteins are 55% and 70%, respectively, and the identity of the active site residues (Ser-Glu-His) is conserved (Fig. S1). Knowledge of the structures and conserved residues of the two enzymes is essential for compound selection and determination of binding modes of the cholinesterase inhibitors.
To validate the robustness of Glide used for the docking, re-docking of the co-crystallized ligands was performed. The crystal structure of AChE (PDB ID: 4M0E 22 ) with Dihydrotanshinone I as ligand and BChE (PDB ID: 5DYW 23 ) with the ligand N-{[(3S)-1-benzylpiperidin-3-yl]methyl}-N-(2-methoxyethyl)naphthalene-2-sulfonamide were selected for redocking. Analyzing the root mean square deviation (RMSD) results for the compounds in AChE and BChE, respectively, gave the values 1.276 Å and 1.570 Å for the two ligands, thus confirming the validity of the methodology used. Moreover, the docked poses of both ligands overlay perfectly with their crystal structure conformations (Fig. S2).
The active sites of holoenzymes of AChE and BChE were defined based on the co-crystallized ligands (Fig. 2B,C) 40 followed by Induced Fit Docking (IFD) of the ligand set containing Ondansetron and the two FDA approved cholinesterase inhibitors Tacrine and Rivastigmine (Table 1). Table 2 shows the docking results and free energies of binding calculated using MMGBSA for both enzymes. The trends between the docking scores and MMGBSA energies are highly consistent, and indicate that Tacrine is the strongest binder in both systems, and Rivastigmine the weakest. Ondansetron seems to bind better to BChE, whereas Rivastigmine is more potent towards AChE. Figure 3A shows the interactions of Ondansetron docked in the active site of AChE. Ondansetron forms π-π stacking interactions with the aromatic rings of Phe295, Phe338, and Trp86. In addition, Trp86 displays one π-cation interaction. Tyr124 and Tyr337 forms hydrogen bonding and π-cation interactions, respectively. There are no direct interactions between the ligand and residues of the catalytic site; however, as outlined above, some aromatic residues in the peripheral site located at the entry to the active gorge are responsible for binding many inhibitors and also contribute to the catalytic efficiency and activity of the enzymes. Among these, Trp86 and Tyr337 are the two major aromatic residues that are involved in binding to the ligands 22 www.nature.com/scientificreports/ Non-competitive binding modes were also observed for Tacrine and Rivastigmine (Fig. 3C). Tacrine forms a hydrogen bond with Asp74, π-π stacking interaction with Phe338, and two π cation bonds with Trp86 and Tyr337. In addition, π-π stacking interactions with Phe338 and Phe295 were observed for Rivastigmine, as well as two hydrogen bonds with Gly121 and Gly122.
In BChE (Fig. 3B), one hydrogen bond and one π-cation interaction are predicted to be formed between Ondansetron and one of the main catalytic residues, His438. The most significant residues to have interactions with ligands in the peripheral site of BChE are Trp82, Trp231, and Phe329 41 . This is in accord with the current study, where we see that Trp82 forms π-cation interaction with Ondansetron, and this amino acid plus Trp231 and Phe329 form π-π stacking interactions with the aromatic rings of the inhibitor. The interaction with His438 in the binding pocket and the residues at the peripheral site would suggest that Ondansetron binds to BChE through mixed inhibition, i.e., binding the active site of the enzyme (competitive inhibition) and the enzyme-substrate complex (noncompetitive inhibition) with different affinity. In comparison, Rivastigmine forms one π-π stacking interaction with Trp82, and three hydrogen bonds with Asp70, Gly116 and Gly117. For Tacrine, there is π-π stacking interaction with Trp82, an ionic bond with Asp70, and a hydrogen bond with the main catalytic residue of His438 (Fig. 3D).
From the free energy of binding values for Ondansetron and the two FDA-approved drugs in AChE and BChE (Table 2), we note that Ondansetron and Tacrine bind most strongly to the holo form of the enzyme. However, to gain better insight on the possible activities of the compounds and the stabilities of the formed complexes, PASS analysis, MD and BPMD simulations were conducted as follows.
PASS prediction. The prediction of activity spectra for substances (PASS prediction) was conducted for Ondansetron with respect to four relevant biological activities. The PASS prediction results of the selected ligand along with the positive controls Tacrine and Rivastigmine are listed in Table 3.
There are three activity levels indicated by the PASS program. If Pa > 0.7, the substance is very likely to exhibit experimental activity; if 0.5 < Pa < 0.7, the substance is likely to exhibit the activity but with less probability; If Pa < 0.5, the substance is unlikely to exhibit. However, if activity is confirmed experimentally despite low Pa value, the substance may be a 'new chemical entity' not catered for in the training of the PASS algorithm 43,44 . As can be seen from Table 3, none of the compounds is predicted to be active, or even moderately active, despite two of them being clinically proven and approved for the treatment of AD. No activity was predicted by this program for Ondansetron as Acetylcholinesterase inhibitor and a very low value was found for Butyrylcholinesterasesimilar to that of Rivastigmine. However, given the experimental evidence (see below), Ondansetron can in this context be labeled a "new chemical entity". A closer look at the table indicates that Ondansetron displays higher active possibility for the treatment of AD in comparison with Tacrine and Rivastigmine. For neurodegenerative disease treatment the PASS data for Tacrine and Rivastigmine are similar (∼ 0.4), and slightly higher than the value for Ondansetron (∼ 0.3). Taken together, it hence seems plausible that Ondansetron is a good candidate for dementia treatment.
Molecular dynamics simulations. The stability and accuracy of the different protein-ligand complexes obtained from the docking studies were explored through 200 ns MD simulations.
The RMSD data for the ligands binding to AChE (Fig. 4A) shows a continuous increase throughout the simulation for Rivastigmine, indicating that this ligand may be more prone to dissociate. We find that, despite the ligand remaining largely in the binding cavity of the protein and still has interactions with some crucial amino acids such as Trp86 in the last snapshot of the MD simulation, the changes are significant with part of the aromatic ring and the alkylamine substituent (Table 1) exposed to the solvent and undergo large fluctuations. This is in agreement with the lower docking score and MMGBSA values ( Table 2) for Rivastigmine. Our results suggest that there are no considerable changes in RMSD values for Ondansetron or Tacrine, and that both form stable complexes with AChE.
The RMSD graphs for BChE (Fig. 4B) show that all three inhibitors have preserved their binding affinity and are still tightly coupled to their respective binding sites. Ondansetron lies below ∼ 2 Å and no significant fluctuations were observed throughout the simulation period, implying that the binding of the ligand to the active site of BChE is very stable and strong. As in AChE, Tacrine in BChE displays slightly lower RMSD values www.nature.com/scientificreports/ compared to our target ligand; however, analysis of the MD trajectory shows that part of the ligand is exposed to the solvent which accounts for the minor fluctuations seen in the graph during the simulation. Analysis of the MD trajectory for Rivastigmine shows that the ligand remained stable in the binding pocket after an initial     We also analyzed the protein root-mean-square-fluctuations (RMSF), showing which residues that fluctuate the most during the MD simulation (Fig. 4C,D). The RMSF values of AChE with any of the ligands bound remain relatively stable (Fig. 4C). The highest variation is related to residues 482-486, when the protein had Ondansetron as ligand. These are part of a loop close to the C-terminal and are not involved with the active site of the protein. The red arrows in the graph indicate the location of the key amino acids in the peripheral site (1, 2, 3, 4, 5 for Trp86, Tyr124, Phe295, Tyr337 and Phe338, respectively). For AChE with Rivastigmine bound, some fluctuations are observed for residues 281-284, and very close to Phe295 (arrow number 3), indicating that the interaction of the ligand with these parts creates changes to the protein. Except for these parts, the protein was highly stable throughout the simulations, irrespective of ligand. We also note that the catalytic residues (black arrows 6, 7, and 8) display very low RMSF values during all the simulations.
The RMSF values for BChE with Ondansetron bound (Fig. 4D) indicate that this complex is quite stable during the whole simulation time, while for Rivastigmine and Tacrine some fluctuations are observed. There are two regions with large peaks close to amino acids in the binding site: residues 70-77 corresponding to a loop including Asp70 in peripheral site (arrow 4) and Trp82 (arrow 5); and residues 334-341 corresponding to a loop very close to Phe329 (arrow 7). For the BChE-Rivastigmine complex, residues 101-105 show large fluctuations; this is a region that is far from the binding pocket. For the system with Tacrine bound, two large peaks are seen near the C-and N-termini, far from the active site; as is also the region with residues 281-285. The location of the three main residues in the catalytic pocket, Ser198, Glu325 and His438, respectively, are shown by black arrows. These residues as well as Trp231 (arrow 6) were found to remain highly stable during the MD simulations of BChE.
Further examination of the system stability was conducted by probing the evolution of interactions between ligand and protein, here measured by the formation of hydrogen bonds during the MD simulations. Ondansetron bound to AChE maintained an average of 1.90 hydrogen bonds, whereas Rivastigmine and Tacrine had an average of 0.45 and 2.47 hydrogen bonds, respectively (Fig. 4E). This indicates that our targeted ligand is more stably bound than Rivastigmine, in agreement with the RMSD plots in Fig. 4A.
In the case of ligands bound to BChE, the average number of hydrogen bonds for Ondansetron is 1.22 while for Rivastigmine and Tacrine it is 1.28 and 2.38, respectively (Fig. 4F). This confirms that Tacrine interacts stronger than the other ligands in this system; for Ondansetron and Rivastigmine the stability and number of hydrogen bonds are very similar.
Binding pose metadynamics (BPMD). Next, we performed binding pose metadynamics (BPMD) simulations for all six systems. The pose stability was evaluated based on the PoseScore, i.e., the RMSD of the ligand with respect to the initial ligand heavy atoms coordinates. The threshold value for ligand stability in the binding pocket is a PoseScore ≤ 2 Å 34 . Figure 5 shows that when Ondansetron is located in AChE the calculated values of PoseScore is 1.10 Å, indicative of stable binding. For Tacrine in AChE, the obtained PoseScore is 1.03 Å, and for Rivastigmine 2.15 Å with a large increase in value, indicating lower ligand stability in comparison with the other two ligands. When Ondansetron is bound to BChE, the PoseScore value is 1.29 Å and for Tacrine and Rivastigmine it is 1.13 Å and 2.39 Å, respectively. Based on these results and the RMSD data (Fig. 4A,B) and the energetics reported in Table 2, we conclude that Ondansetron can be a more potent candidate against cholinesterase compared to Rivastigmine; however, less so than Tacrine in both proteins. BChE activity was investigated through different in vitro assays, and the type of inhibition and the kinetic factors of the binding calculated. The enzymatic activities of AChE and BChE were evaluated spectrophotometrically at room temperature using Ellman's method as described in the method section, with the rate of increase in absorbance at 410 nm followed for 5 min. The experiments were done in triplicate. The tested concentrations of AChE and BChE ranged from 0-84 µM and 0-42 µM respectively. Lineweaver-Burk plots were used to determine the type of inhibition. It was found that the drug inhibits AChE by non-competitive inhibition (Fig. 6A), while for BChE Ondansetron exhibits mixed inhibition (Fig. 6B). The data thus confirm the analysis from the docking regarding the type of inhibition (Fig. 3). The K m of AChE was constant in the presence and absence of the drug and was calculated to ∼ 0.27 mM, while the K m of BChE was variable depending on concentrations of the drug.
Affinity of Ondansetron in the enzymes. Since the K m of AChE in the presence of the drug was constant, V max was used to calculate IC 50 in this case (Fig. 7A). For BChE, K m of mixed inhibition was instead used to determine the IC 50 value (Fig. 7B) Ondansetron is hence a more potent binder to AChE than Rivastigmine, but less so than Tacrine, which agrees perfectly with the MMGBSA data of  46,47,48,49 , indicating that Ondansetron is at least 3-8 times more potent than Rivastigmine but again less effective than Tacrine. Again, this trend is fully supported by the computed binding energies (Table 2). We note however that the most important adverse effect of Tacrine is liver toxicity, gastrointestinal problem, nausea, and vomiting and that it due to this has been withdrawn in some countries 11 . We thus propose that Ondansetron may be a more potent alternative drug targeting cholinesterase than Rivastigmine, and with relatively mild side effects as compared to Tacrine. Calculation of the inhibition constants (K i ) were carried out through Lineweaver-Burk and their secondary plots (Fig. 7C,D). The x-intercept indicates the K i value, giving K i = 35 µM and 6.1 µM for Ondansetron in AChE and BChE, respectively, which also confirms the higher binding affinity of Ondansetron in BChE noted computationally.
Effect of pH on Ondansetron binding. Both cholinesterases are pH dependent and perform optimally at pH above 7 50 . In general, the rate of enzymatic hydrolysis of acetylcholine and butyrylcholine increases from pH 6 to 9 and shows maximum activity at pH 8. However, the activity decreases markedly at pH 11 51 . Several physiological and pathophysiological conditions contribute to the pH in the neurons. Both an ionic site and positive groups on the surface of ACh receptors recognize this protein for the attachment (Fig. 1) as a result both groups are involved in modulating the pH value in neuron cells. Moreover, some selective ions in the membranes of nerve cells contribute to fixing alkaline charges 52 .
To this end, the activity of AChE and BChE was investigated within the pH range 5-11 at 30 °C, using 21 µM Ondansetron. The optimum activity for AChE alone was observed at pH 9, while in the presence of the drug, a shift to pH 10 was observed (Fig. 8A). For BChE the optimum pH is 9, which remained unaltered in the presence of the drug (Fig. 8B). The pH profiles hence suggest that both enzymes are more active in alkaline pH than in acidic pH, and also show a clear impact of drug binding on the enzymatic activity.
Effect of temperature on Ondansetron binding. The effect of different temperatures (0-70 °C) on the enzymatic activity was also measured in the presence (21 µM concentration) and absence of Ondansetron. AChE showed maximum activity at 40 °C and 50 °C in the absence and presence of the drug, respectively (Fig. 9A). BChE showed the same pattern as AChE, with optimum activity at 40 °C while in the presence of Ondansetron it shifted to 50 °C (Fig. 9B). Both enzymes were also inactivated at 70 °C both in presence and absence of the drug. www.nature.com/scientificreports/ Arrhenius plots were obtained to calculate the activation energies of the enzymes and the effect of the drugs thereon. It indicated that, in the presence of Ondansetron, more energy was needed for the reaction to proceed, suggesting that binding of the drug to the enzymes will increase the activation energy. As a result, Vmax of the enzymes were calculated at different temperatures in order to generate Arrhenius plots. For AChE, Vmax increased from 25 to 30 kJ mol −1 in the presence of the drug (Fig. 9C) while for BChE the activation energy was 23 and 32 kJ mol −1 in the absence and presence of the drug respectively (Fig. 9D). Again, this supports the finding of Ondansetron binding stronger to BChE, as noted in Table 2.

Conclusions
Ondansetron has herein for the first time been investigated as a potent inhibitor for cholinesterase proteins. In the current study, computational modelling and experimental studies were carried out to explore structural properties and binding mechanisms of the compound. Our observations suggest that it can inhibit the enzyme  www.nature.com/scientificreports/ activity by interaction with residues near or in the active site of the enzymes, in non-competitive or mixed inhibition manners for AChE and BChE, respectively. Analyses of MD simulations revealed that Ondansetron is a stable ligand in both enzymes, and binds with higher affinity in comparison with Rivastigmine. We also carried out different experimental studies by calculating the type of inhibition and K i values from Lineweaver-Burk plots that fully confirm the data from the in silico modeling. The effect of binding to AChE and BChE were investigated, that show that our target has stronger inhibitory activity towards BChE than AChE (IC 50 : 2.5 µM and 33 µM, respectively). However, in relation to the two FDA approved drugs Rivastigmine and Tacrine, and considering the serious side effects such as hepatotoxicity for Tacrine, we herein propose that Ondansetron may be a potential repurposed candidate drug to reduce or block cholinesterase activity. Further experimental studies (in vitro and in vivo) would be required to verify the pharmacological relevance of the observed cholinesterase inhibition by Ondansetron.